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Non-dispersive wave packet for massless fermions is formulated on the basis of squeezed coherent 
states that are put in a form of common eigenfunction for the Hamiltonian and the helicity operator, 
starting from the Dirac equation. The wave packet thus constructed is demonstrated to propagate 
at a constant velocity as that of light. This explicit expression of wave packet for the massless 
fermions can facilitate theoretical analysis of problems where a wave packet is of formal significance. 
Furthermore, extensive wave packet may result in a superluminal velocity statistics if determined 
from the time-of- flight measurement, as recently done on muon neutrinos, when a threshold particle 
flux or energy transfer, which is eventually referred to the propagation of wave packet, to invoke a 
detection event is assumed. 
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In fundamental quantum theory, a free particle is rep- 
resented by a plane wave e 4 ( fca: ~ w *) which is the eigen- 
function of the momentum operator. However, the plane 
wave cannot be normalized, thus it does not fit our pic- 
ture of particles which at least should make their presence 
only in a limited region of microscopic scale Packet 
waves have been formulated for light coming in indivisible 
unit of photons, a kind of massless boson, and for massive 
fermions such as electron, based on a dispersion relation 
derived from the electromagnetic wave equation for the 
former and from the Schrodinger equation or Dirac equa- 
tion for the latter, respectively. Remarkably, however, 
there is no mathematical formulation of wave packet for 
massless fermions, to the best knowledge of the authors. 
With the absence of massless fermions, the picture of 
particles as a wave packet is obviously incomplete. 

With the introduction of wave packet, various particle 
behaviors are often discussed as a matter of fact in terms 
of this concept. Yet in many circumstances, for instance 
when talking about the collision or detection by a detec- 
tor, the particles are simply thought of as a point object. 
The handling of particles as a point object is a valid con- 
venience when the size of the particles is much smaller 
than the scale of the experimental setup, or to be nearer 
to the reality,less than the size of atoms or atomic clusters 
that constitute the responsive agent of a detector. This 
is not a challenging requirement for electrons or protons 
since these massive fermions have a very short Compton 
wavelength and they just behave as distinct identities 
though formally they can also be represented as a wave 
packet. For such massive fermions, the concepts of point 
particle and wave packet do not cause much confusion, 
and the output of measurement for one single particle 
striking a macroscopic detector should give a sharp value 
within the precision of the detection method. In the case 
of massless bosons, things become complicated. Though 



we know that light comes in indivisible unit of photons, 
our description of relationship between light and photon 
is incomplete since, for example, the light can be po- 
larized but we don't know how each photon can keep a 
record of the light's polarity. This would take us away 
from the traditional ideal of elementary particles, and 
more significance may be attributed to the concept of 
wave pack. Put this fact aside, we notice that the pho- 
tons now can be brought into light pulses of enormously 
differing durations. It is reasonably conceivable that the 
wave packets describing light can in principle assume a 
macroscopic extension, noting that a temporal duration 
of 1.0 ns corresponds to a spatial extension of 30 cm for 
a light pulse subject to none of the decoherence mecha- 
nisms. 

A quite natural question may be posed about mass- 
less fermions: how could they be represented as a packet 
wave, and what the extensive wave packet may imply 
for the investigation of behaviors of such particles? This 
may be a question of some relevance, say, with regard to 
the recent report of superluminal velocity measurement 
on muon neutrinos [2j. Neutrino is a massless fermion 
in the standard model. Some experimental observations, 
e.g., neutrino oscillation, point to the probability that 
neutrinos may have mass. It is believed that the upper 
limit for heavier neutrinos is at the order of magnitude 
of 2 eV, and the lightest neutrino could in principle be 
massless 0]. We shall always presuppose m„ = with- 
out distinguish the flavors of the neutrinos in the cur- 
rent work. Now, a massless neutrino, following relativis- 
tic quantum mechanics, should propagate at the velocity 
of light in vacuum, c, and the detection of its velocity, 
if the neutrino behaves as a point particle, is expected 
to give a sharp value of c within the precision of the 
method. However, Adam and coworkers recently showed 
that the velocity of muon neutrino with the average en- 
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ergy of 17 GeV, by time-of-flight measurement, displays 
surprisingly a skewed distribution biased at the > c side. 
This has been even interpreted as the proof of supcrlu- 
minal existence. A detailed inspection of the data for the 
experiment performed at the average of 28.1 GeV further 
indicates that the data taking only events at energy >20 
GeV give to a more severe superluminal velocity statis- 
tics [2j. Could it be possible that the massless neutrinos 
propagate as a train of wave at a group velocity of c, yet 
the behavior of the wave packet could lead to an energy- 
dependent velocity statistics deviating from c? 

In this article, we report the formulation of wave 
packet, starting from the Dirac equation, for massless 
fermions on the basis of the squeezed coherent states [1] , 
and the wave packet thus constructed is non-dispersive 
and propagates at the velocity of light. We see that 
by assuming a threshold particle flux or energy transfer, 
which to the end is referred to the probability distribu- 
tion described by the wave packet, is desired to invoke 
a detection event for the particle, the measured velocity 
statistics for neutrinos could find a reasonable explana- 
tion without breaking relativity. 

For a massless fermion, the Hamiltonian is given in the 
form 



H 



(1) 



Where, c, the velocity of light, is a parameter that con- 
nects space and time. This doesn't necessarily imply that 
the particle has a definite velocity as that of light, bearing 
in mind the Dirac equation for the massive electron. 

For simplicity, without loss of generality, we con- 
sider the situation that the particle propagates along x- 
direction, then the Hamitonian (1) is rewritten as 



H = ca x p x 



(2) 



Now, let's represent the position operator x and momen- 
tum p in terms of the creation and annihilation operators 



tf W _ + 4) . 



(3) 



With A being a parameter in the dimension of length, 
then the Hamitonian in Eq.(2) turns into 



H = Aa x (b + tf), 



(4) 



Where A = 



V2A ' 



It is obvious that the helicity, r\ = ~& ■ ~P /p, i.e., the 
projection of the spin of the particle onto the direction of 
momentum, is a conserved quantity since it commutates 
with the Hamiltonian in Eq.(2). The eigenstates for the 
helicity operator are given as follows 

_ / (e-^+^IO) \ 
\+> y (e -|6t 6 t +fe t )|0) J , 

_ ( (e-i» tb '- bt )\0) \ 



(5) 



which correspond to the eigenvalue +1 and -1, respec- 
tively. The eigenstates in Eq.(5) are, however, not 
normalizable since the integral arising from the term 



/3f> T f> T 



|0) is unlimited for 
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Consequently, such 



eigenstates of a unique helicity is not proper for the phys- 
ical representation of a massless fermon though it was 
claimed that the state of a massless neutrino could never 
be a superposition of two different hclicitics and the ex- 
periment gave a direct indication that the helicity of neu- 
trino is negative [f| . In order to construct a normalizable 
physical state in the form of wave packet for a massless 
fermion, the mixing of helicity seems unavoidable. Let's 
take the wave packet given below as the initial state of 
the massless fermion 



\t = 0) 
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da, (6) 



where the weight factor in the integral is Gaus- 
sian. It can be shown that such an initial state 
is normalizable, as can be proved that (t = 0\t = 

A 2-y 2 

0) = ^7y= 1 e 4 --i- 1 °' , in deriving this we have used 
the enclosure relation for the coherent state |$) = 
e^ + |0), i.e., / d(Re£l4(Im£l e _(Re<E. 2 +Im$ 2 )| $ ^ $ | = 

1, and the integrals e~^ x2 s\a(bx)Ax — and 



X^oo e 13x2 cos (bx)dx — . /¥ 



e 4 * 3 . 



The temporal evolution of the state represented by the 
wave packet in Eq. (6) can be derived in a straightforward 
way. At any moment t, it has 
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e -j{a-a ) 2 e -iE{a)t 
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da. (7) 



Thus the expectation value for the position operator 
at any moment t is 



x(t) 



_ {t\f 2 {b-tf)\t) 

(t = 0\t = 0}~ (t = 0\t = 0) : 



(t\x\t) 



(8) 



which leads to the result x(t) = V2A • At, implying that 
the velocity of the packet wave is a constant of c, as 
t, = m. = V2A • A = c. 

Now let's calculate the width of the wave packet at 
any moment t, V = {(x{t) 2 ) - (i(i)) 2 ) 1 / 2 . Note that 



-4-(t|(6 - 6+) 2 |t), it thus has (x(t) 2 ) 



_^J±_ = 2A 2 (A 2 i 2 + I) - f , consequently 

r(t) = [2A 2 (A 2 t 2 + -)- — -(V2AAt) 2 } 1 / 2 = A{---) 1 ' 2 . 

7 2 7 2 

(9) 

Based on the arguments above one concludes that the 
wave packet given in Eq.(6) is non-dispersive, and prop- 
agates at a velocity of c, therefore it can be used to repre- 
sent a massless fermion or a beam of massless fermions.lt 
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is noteworthy to mention that 4 > 7 > 1/4 as required 
by the normalization condition and the positiveness of 
the wave packet width. 

From Eq.(9) we see that the extension of this non- 
dispersive wave packet is proportional to the free param- 
eters A, consequently it can in principle be made suffi- 
ciently extensive. This is analogous to the case of light for 
which the wave packet can be of an enormously variable 
lengths, since the Compton wavelength of photon is in- 
finitely long and the photons are almost decoupled from 
each other. For massless particles such as photons and 
neutrinos, their behaviors are expected to display more 
undulatory characteristics. In contrast, massive particles 
such as electrons and neutrons have a limited Compton 
wavelength, their wave packet can have only a limited 
extension-in general cases they are just treated as a dis- 
tinct particle instead of a train of wave. 

Now it is probably appropriate to discuss the statis- 
tics of velocity measurement for the massless fermions, 
e.g., for neutrino. When the neutrinos are represented 
by the wave packet formulated in Eq.(6), does that imply 
the possibility of measuring a superluminal velocity un- 
der certain circumstances? Or, what kind of an energy- 
dependent velocity measurement can be anticipated for 
neutrinos when taken as a non-dispersive wave packet 
with a constant propagation velocity of c? 

The difference of the measured velocities for muon neu- 
trinos relative to that of light has been given in a statis- 
tical notion, and the original data of velocity scatter ran- 
domly around the velocity of light. For the muon neutri- 
nos with an average energy of 17 GeV, Adam and cowork- 
ers, based on the measurement of time-of-flight through 
a base line of 730 km long, came to the conclusion that 
(v-c)/c= (2.48±0.28(stat.)±0.30(sys.))x lO- 5 0. A 
measurement performed for an average neutrino energy 
of 28.1 GeV shows only negligible difference from the 
result above, but the statistics taking events of energy 
higher than 20 GeV shows a significantly larger velocity 
than that for energy below 20 GeV, the early neutrino 
arrival time St is ~67.1 ns in the former case while only 
~53.1ns in the latter case Q. Furthermore, it is worth 
noting that the Minos experiment with neutrino energy 
peaking at ^3.0 GeV with a tail extending to above 100 
GeV results in (v - c)/c = 5.1 ± 2.9 x 10~ 5 >], but a 
stringent limit of (v - c)/c < 2 x 10~ 9 (given in abso- 
lute value!) was set by the observation of (anti) neutri- 
nos emitted by the SN1987A supernova in the 100 MeV 
range @, S]- Roughly speaking, the measured neutrino 
velocity seems to suggest that following a larger neutrino 
energy there may be an increased probability of measur- 
ing a superluminal velocity. 

The wave packet formulated in Eq.(6) for massless 
fermions is strictly non-dispersive, and it propagates at a 
constant velocity of light. This concludes from the solu- 
tion of the equation of motion for the particles, which is 
consistent with both special relativity and quantum the- 



ory. If the velocity of such a wave packet is determined by 
measuring the time of flight, then so long as time tagging 
is performed with respect to precisely the same position 
on the wave packet, say the center of the wave packet but 
not necessarily, it will result in the same value, c, even 
though it has a finite extension. For the muon neutrinos 
in Ref.Q, it requires a time tagging referring to exactly 
the same position of the wave packet at the source and 
at the detector, a quite challenging task because the time 
tagging was performed via different processes. The finite 
extension of the wave packet may cause a peculiar sys- 
tematic deviation of measured neutrino velocity from the 
velocity of light, following the Copenhagen interpretation 
of the wave function. 

The event that a particle initiates as being registered, 
whether on the surface of a solid detector directly or via 
an intermediate nuclear reaction, is essentially quantum 
mechanical, i.e., in a statistical sense. Suppose that it 
needs a threshold particle flux, or energy transfer, to in- 
voke a detection event, e.g., to generate taus to be mea- 
sured to indicate the arrival of the neutrinos Q , then this 
implies that a threshold portion of the wave packet has 
to enter the detection region or the device @, [l(| ■ As we 
mentioned above, the wave packet for neutrinos can have 
a macroscopic extension as a pulse of light, the super- 
luminal velocity statistics for neutrino can find a natu- 
ral and straightforward explanation under the threshold 
probability assumption. In case a larger threshold por- 
tion of the wave packet, implying a larger particle flux or 
energy transfer, is required to invoke a detection event 
that the time tagging at the detector might be retarded 
relative to the time tagging at the source, the velocity dis- 
tribution will be biased by < c results; and if a smaller 
threshold probability of the wave packet is required that 
the time-of flight measurement may lead to an early ar- 
rival time, then the velocity distribution may be biased 
by > c results (Figure 1). Clearly, for the beam of neu- 
trinos of greater energy, the velocity distribution may 
contain more superluminal data, farther from the value 
of c. 

In summary, non-dispersive wave packet with a con- 
stant propagation velocity of c has been formulated for 
massless fermions starting from the Dirac equation. Such 
a wave packet is constructed from the superposition of 
common eigenstates of the Hamiltonian and the helic- 
ity operator, which are represented on the basis of the 
squeezed-coherent states. Based on the assumption that 
a threshold portion of the wave packet is demanded to in- 
voke a detection event, the wave packet, macroscopically 
extensive as it may be in principle as a beam of light, 
can provide a quite reasonable explanation to the super- 
luminal velocity statistics recently reported for neutrino, 
which is anticipated to occur more frequently and farther 
away from the velocity of light at higher particle energies. 
The result can serve as a basis for the discussion of other 
problems concerning massless fermions where the formu- 
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FIG. 1: Velocity measurement of a non-dispersive wave packet 
for massless fermion with a constant propagation velocity of c. 
In case that the portion of the incident wave packet demanded 
to invoke a detection event is (a) ahead of, (b) coincident with, 
or (c) lagging behind the point registered as t = , then the 
velocity calculated from the time of flight will be v > c, v = c 
or v < c, respectively. 



lation of the wave packet may be of formal significance. 
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APPENDIX 

1. The module square of the initial state \t = 0) 
For 



0) 



-7(0— q )- 
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da, 



(t = 0\t = 0) = / e- 7(Q - Qo)2 - 7(a ' 2 ~ Qo)2 dada' • 2(0|e^ bfc+Q V5 fctbt+Q ' bt |0) 



= 2 J e -^ a - ao)2 -^ a ' 2 - ao)2 dada J d ( Re$ ) d ( Im$ ) c -2Re$ 2 +faW)Re$+»f Q ~a , )Im$ 
For brevity, replacing Re$ with £ and Im$ with r), 

(t = 0\t = 0) = - y" e -7(^"o) 2 -7(a' 2 -a ) 2 dada ' y" d ^ d7?e -2C 2 + (a+a')C+*(a-a')r ) 
2 
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e -7(a-ao) 2 -7(a' 2 -a ) 2 dQ;da ' / d £ d?7e -2C 2 + (« W K [ cos ( a _ q,')^ + j sin ( a _ q/)^]. 



Discarding the imaginary part of the integral as it must vanish, one obtains 

(t = 0|t = 0) = ^ ye- 7(a - Qo)2 - 7(Q ' 2 - Qo)2 dada'ydCdr/e" 2c2+(a+Q ' )c cos(a-a')77 

= \ J e-^ a ~ ao)2 -^ ia ' 2 - ao)2 dada J dCdr]e- 2i ^^ )2+i ^^ cos{a-a )r) 
= - I e -7(a-«o) 2 -7("' 2 -"o) 2 ! dada ' I dC'd?7e- 2C ' 2+i2± ^ Ji cos(a-a')7 ? 



~V2 / e ~ 7(Q ~ Qo)2 ~ 7(Q ' 2 ~ Q ° )2dadQ/ / drje 1 ^ 1 cos(a-a')rj 
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= \[\ J e -~ ,{a2 - 2aaa+a « )+ ^ + ^ e~^ a ' '-^ f+^ COS (a - a )r)&a&d dTj 
= \[^z~ ial y'e-^-5)« 2 +^^« e -7(a'-«o) 2 +^ cos ( a _ Q ,' ) r?dQ , dQ ,' dr7 

^e-^o f e -h-i^" 2 cos(a " + 8 ^"+ l Q - a ' J^e^-^'-^^^da' d?7 
n J 8(7 -i) 



2 8 7 2 -2 7 2 r i. "2 // 870:0 - (87 - 2)a 

-e s -r-i ai e n » >a cos a 7700s— 1 — — 77 

7r 7 87 - 1 

sina^sin 8700 zfel 2 )" ^^H^-^)*' W«o«'W ' da d?? 



2 87^27 2 /• / 8tt 2_ n 2 8700 - (87 - 2)a 

-e 8 7 -i a ^ / ,/- -e 87-1^ cos— 1 „ — 7? 

y v 87 - 1 87-1 ' 

«7 2 -27 f ' 87 2 <»0 x2 M ^°| ) , 

X e 87-1 8-, , ->-2-,> e (87-l)(87 2 -27 ; dQ; d? 7 

16 4 ,7 2 2 /" 2 „2 s T 2 -2 T '"2 . 870-0 87— 2, /// 87 2 o Nn 

e i7^^ Q o / e ~— v e ~^^r<* cos[ — - — a + — ^ — — ) ryda d?7 

87-I 7 '87-I 87-I 1 87 2 - 2 7 /J/ ' 
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2. Motion of the wave packet 

For state \t = 0), its temporal evolution is given by 



!*> = /■ 



e -i6t b t +a6 t |0) Ida- Uij 



Thus it has • (t = 0\t = 0) = (t\^(b - tf)\t) 
from which, following the same tricks in deriving (10), 
one obtains 

x ■ (t = 0\t = 0) = V2AM - ^ e^ a °. (12) 

V47- 1 

This implies that x(t) = V2AAt, consequently the veloc- 
ity of the wave packet is 



v=^- = V2AA = c. 
3. Dispersion of the wave packet 



(13) 



The width of the wave packet is given by T = 
[(t\(x - x) 2 |/j)/(7j||7j)] 1 / 2 which can be simplified to T = 
_»L. _ a^i/^ Since 



L(t=0|t=0) 



iA 



(t\x 2 \t) = (t\(—(b-tf)f\t) 

= + fo t2 - 26t6)|t> -¥.( t \tHU) 

Thus we just need calculate (t\x 2 \t) + ^-{t\t) which, 
again following the similar techniques in deriving (23), 
can be found to result in 



(t\x 2 \t) + ^-(t\t) 



8A 2 



1 



_, r - e ^^ a °(A 2 t 2 + -), (15) 
V47 - 1 7 



or = 2A 2 (A 2 t 2 + I) - w hich further leads to 

T(t) = A(- — I) 1 / 2 , i.e., the width remains constant. 
The wave packet is non-dispersive. 



